Dynamics of a domain wall in a magnetic nanostrip: a toy model 
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In this report we demonstrate a simple model for the motion of a vortex domain wall in a fer- 
romagnetic strip of submicron width under the influence of an external magnetic field. The model 
exhibits three distinct dynamical regimes. In a viscous regime at low fields the wall moves rigidly 
with a velocity proportional to the field. Above a critical field the motion becomes underdamped 
as the vortex moves periodically across the strip; these oscillations are accompanied by a slow drift 
with a decreasing velocity. At still higher fields the drift velocity starts rising linearly with the field 
again but with a much lower mobility dv/dH than in the low-field regime. We calculate the relevant 
quantities and compare them to experimentally observed values. 



I. INTRODUCTION 

Dynamics of domain walls in ferromagnetic strips and 
rings with submicron dimensions is a subject of active 
researchjii^i^ This topic is directly relevant to several pro- 
posed schemes of magnetic memory and is also interesting 
from the standpoint of basic physics. The dynamics of 
domain walls under an applied magnetic field has distinct 
regimes: viscous motion with a relatively high mobility 
at low fields and underdamped oscillations with a slow 
drift at higher fields.'' 

The nontrivial dynamics is related to the composite 
nature of a domain wall in a nanostrip: it consists of a 
few — typically two or three — elementary topological de- 
fects in the bulk and at the edge of the strip. As a re- 
sult, a domain wall has several low-energy degrees of free- 
dom that are relevant to the dynamics. Weak external 
perturbations engage only the softest (zero) mode — rigid 
translations along the strip. Larger external forces excite 
higher modes thereby altering the character of motion. 

The general approach to the dynamics of domain walls 
in thin ferromagnetic strips with a submicron width w 
and thickness t ^ w was described recently by Tretiakov 
et al.^ The configuration of a domain wall is parametrized 
by a few collective coordinates ^ = {£,i,£,2, ■ ■ ■ ,£,n} and 
the free energy of the system U is treated as a function 
of ^. The resulting equations of motion can be written 
in the vector notation as 

F - + = 0. (1) 

Here components of the vector F are generalized forces 
Fi = —dU/dS^i, the symmetric matrix F and antisymmet- 
ric matrix G represent the viscous and gyrotropic tensors, 
respectively. 

The main goal of this paper is to illustrate the 
collective-coordinate approach on a very simple model 
of a vortex domain wall^ that served as a prototype for a 
more realistic model of Youk et alJ Despite its simplic- 
ity, the model captures all of the main features of a vor- 
tex domain wall and yields simple analytical results for 
the relevant physical quantities. Quantitatively speaking, 
the values of the forces computed in this model deviate 
by no more than 30% from those obtained in the more 



realistic model of the vortex wall. Thus one can make 
meaningful comparisons between the analytical results 
obtained in this paper and experimental data. 

In the main body of the paper we describe the simpli- 
fied model of the wall and compute the generalized forces 
and the viscous and gyrotropic tensors. By substituting 
these quantities into Eq. ([T]) we obtain simple equations 
of motion. At low fields the equations describe steady 
viscous motion of the wall with a velocity proportional 
to the applied field. The vortex is shifted in the trans- 
verse direction by an amount proportional to the velocity 
of the wall. At a critical velocity the vortex is expelled 
from the strip and the steady motion breaks down giv- 
ing way to an oscillatory regime. As the applied field 
increases further, the drift velocity decreases at first but 
then again becomes proportional to the applied field; the 
mobility coefficient dv / dH is substantially lower than the 
corresponding value in the viscous regime at low fields. 
These results are compared to experimental data. 



II. MODEL WALL 

In our calculations of the wall dynamics, we use a sim- 
ple model of the vortex domain wall consisting of four 
domains with uniform magnetization and separated by 
90° Neel walls (Fig. [J). 

We assume that only two softest modes of the vor- 
tex wall arc involved in magnetization dynamics, so that 
the configuration of the wall is fully described by the 
two coordinates {X, Y) of the vortex core. In that case 
the equations of motion ([T]) refiect the balance of forces 
acting on a particle moving in two dimensions with a 
velocity V — {X,Y). The forces include a conserva- 
tive term {—dU/dX,—dU/dY), a viscous term FV, and 
a gyrotropic term^ GV = pGz x V. The gyrotropic 
force depends on the out-of-plane polarization of the vor- 
tex core p = M™''®/|M™''''| and the gyrotropic constant 
G — 2TTJt, where J = fioM/j is the density of angu- 
lar momentum, - The conservative and viscous terms are 
discussed next. 
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FIG. 1: A simple model for a vortex domain wall.^ The pan- 
els show states with a fixed longitudinal coordinate of the 
vortex X = const; the transverse coordinate is Y = w/2, 
w/4, 0, — w/4, and —w/2. The vortex core is denoted by the 
filled circle. Shaded areas indicate the locations of magnetic 
charges. 



A. Free energy U and conservative forces 

In strips that support vortex domain walls, the domi- 
nant contribution to the energy in the absence of an ap- 
plied field is due to magnetostatic interactions. For any 
domain wall in a strip of width w, thickness t, and satu- 
ration magnetization M, there is a total magnetic charge 
Q = 2fioMtw associated with the wall.^ In a vortex wall, 
nearly all of this charge is expelled to the edge. In our 
simple model, magnetic charges form two lines of lengths 
w — 2Y and w + 2Y with constant density of charge per 
unit length p — /igA/i. 

The magnetostatic energy of this wall E{Y) = E{0) + 
+ 0{Y^) has a minimum at F = 0. This leads to 
a force —kY that acts to keep the vortex centered on the 
strip. In general, the total exchange energy of the wall 
may change with the position of the vortex, altering the 
restoring force slightly. However, in our simplified model, 
the exchange cost comes entirely from the four Neel walls 
that make up the vortex wall. Because the length of these 



walls does not change as the vortex moves, we need not 
consider the exchange interaction in our analysis of the 
wall dynamics. 

A line of charges of length L has the self-energy 



Eo{L) 



Stt 



(2) 



The divergence a,t x — x' requires a regularization. In a 
crude way this can be done by introducing a short-range 
cutoff in the integral, i.e. by integrating over distances 
jx — a;'| > Ci, where C is a numerical constant. We then 
obtain a logarithmic dependence on Ct: 



EoiL) 



in 



■L[\og{L/Ct)-l]. 



(3) 



The self-energies of the two lines of charge is 



EselfiY) 



Eoiw - 2Y) + Eoiw + 2Y) 
fioMH'^Y^ 



Escu{0) 
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+ OiY^). (4) 



Note that the cutoff Ct affects only the constant term 
Escu{0)', the quadratic term is not sensitive to the exact 
value of C. 

In a similar way we evaluate the interaction of the two 
lines of charges. 
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The sum of the quadratic terms in Eqs. (|4]) and ^ yields 
the "spring" energy kY"^ /2, from which we determine the 
spring constant: 



k 



2(1 - fioMH^ 



(6) 



Next we deal with the Zeeman energy of the wall 
—Hot J d'^x H • M in the presence of an applied magnetic 
field H parallel to the axis of the strip. A longitudinal 
shift of the vortex by AX results in a pure translation 
of the wall. Independently of the wall shape, the rigid 
shift changes the Zeeman energy by —QHAX, where 
Q = 2p,oMtw is the magnetic charge of the wall. There- 
fore the longitudinal Zeeman force is QH in any model. 

The Zeeman force also has a transverse component. 
As can be seen from Fig. [1] transverse motion of the vor- 
tex core changes the total magnetization AI^ of the strip 
and thus affects its Zeeman energy. As the vortex core 
crosses the strip from top to bottom (Fig.[T]), the Zeeman 
energy decreases linearly by ^^.^HMtw^ . Therefore the 
transverse component of the Zeeman force is —2QH. 

The total free energy of a wall with the vortex core at 
{X, Y) is thus 



U{X, Y) = kY^/2 + 2QHY - QHX. 



(7) 
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B. Viscosity tensor f and viscous drag 

We next consider the viscosity of the vortex waU. The 
viscosity that appears in Eq. ([T]) is a symmetric matrix 
whose components are given by)^ 



r. 



aJt d X 



(8) 



where (f> is the azimuthal angle characterizing magneti- 
zation. 

An infinitesimal shift in the collective coordinates X 
and Y affects magnetization in the vicinity of the Neel 
walls only. We begin by considering the contribution 
of a single Neel wall emanating from the vortex core at 
±45°, (t){x,y,X,Y) = f{x-XTy±Y). For such a wall, 
derivatives with respect to collective coordinates can be 
reduced to ordinary gradients: d(j>/dX — ~d(j>/dx — — /' 
and d(p/dY = —d(f)/dy = ±/'. As a result, the tensor 
components are equal to each other, up to a sign: 



XX 



YY 



XY 



aJt / (fxf 



Note that this represents, up to a trivial constant, the 
exchange energy of the Neel wall, which has been calcu- 
lated, e.g., in Ref. [o!. We thus obtain viscosity coeffi- 
cients for the Neel walls intersecting at the vortex core, 
Txx ~ Tyy = T^xy = 0.152aJtw/ \, where the ex- 
change length A — ^ Aj yn^M'^ — 3.8 nm in permalloy. 
Opposite signs of the off-diagonal component Vxy can 
be understood by noting that, as the vortex moves along 
y, the two Neel walls shift along ^X and —X creating 
equal and opposite viscous forces in the X direction. 

The two peripheral Neel walls have the functional form 
(j){x,y,X,Y) = f{x + y — X + Y ± w), so that their 
contributions are the same as that of the central Neel 
wall perpendicular to them. Adding the contributions of 
all four Neel walls yields a total 



Txx = Tyy = -2TxY = 0.608aJiw/A, 



(9) 



independently of the vortex position. 

It is instructive to compute the ratio of the viscous and 
gyrotropic forces: 



XX 



/G = Omiaw/X. 



(10) 



The small value of Gilbert's damping in permalloy, a ~ 
0.008,^" leads to the dominance of the gyrotropic force in 
strips with submicron widths. The smallness of Txx/G 
can be exploited to organize an expansion in powers of 
this small parameter. 



III. WALL DYNAMICS 

Equations of motion ([T]) for two generalized coordi- 
nates = X and £,2 —Y read 



where Fy = Tji is a viscosity tensor, p is the polarization 
of the vortex core, and e^j is the 2x2 antisymmetric ten- 
sor with ei2 — +1^ The generalized forces Fi = —dU/d£,i 
are derived from the free energy ([7]). We thus arrive at 
the following equations of motion for the vortex core: 



X 
Y 



QH , k{TxY-pG) 



Txx dctF 
det F + G2 



(Y-Y,^)^ 



(12) 



where the equilibrium Y position of the vortex is given 

by 



fcFeq = -pG^(l+w), 

i XX 



(13) 



where g = {2Txx + Txy)/G <C 1. It is worth noting 
that the magnitudes of the transverse displacement |l^q| 
are slightly different for the two values of the vortex po- 
larization p. This effect can be traced to the lack of the 
reflection symmetry y i—^ —y in a vortex wall, which leads 
to nonzero transverse components of the Zeeman force 
—2QH and the viscous force TyxX. As a result, trajec- 
tories of vortex cores with p = +\ and —1 are slightly 
different. 

Analysis of the equations of motion yields three dis- 
tinct regimes (Fig. [J). Below a critical field He we find 
steady viscous motion with a high mobility p, = dV/dH. 
Immediately above the critical field He the motion ex- 
hibits an oscillatory component; the drift velocity quickly 
decreases as the applied field grows. At much higher 
fields, H ^ Ho, the drift velocity rises linearly again 
but with a much lower mobility p than at low fields. 
The separation of scales He and Hq is guaranteed by the 
smallness of the parameter Txx /G. 



A. Low field: H < He 

In a low applied field the wall exhibits simple viscous 
motion. The transverse coordinate of the vortex will 
asymptotically approach its equilibrium position Y^q, so 
long as the latter is within the strip. The wall then moves 
rigidly with a steady longitudinal velocity 



X - QH/Txx- 



(14) 



Fi - Fh^j -I- pGe^j^j = 



(11) 



Experimental data of Beach et al^ yield 
Q/Txx ~ 25 (m/s) Oe^^ at low fields for a strip 
600 nm wide, which gives Txx/G = 0.13. while our 
Eq. (HOl) yields Txx/G = 0.12 if we use the value of 
a — 0.008 measured by Freeman et alJ^ 



B. Critical field: H = He 

The low-field regime ends when the equilibrium posi- 
tion of the vortex core is pushed outside the strip edge. 
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FIG. 2: Drift velocity vs. applied field curves for permalloy 
strips of various widths. We use t — 20 nm, a — 0.01, and 
M = 8.6 X 10^ m-i A. 



strip, a new vortex is injected at the same location with 
the opposite polarization —p. The vortex thus moves 
between the edges switching its core polarization each 
time it reaches an edge. 

Once the transverse coordinate of the vortex Y be- 
comes a dynamical variable, the motion acquires an en- 
tirely different character. As we already pointed out, 
the gyrotropic force GV dwarfs the viscous one, TV, in 
permalloy strips. To zeroth order in Txx/G, the dynam- 
ics is conservative: the vortex core moves along equipo- 
tential lines U {X, Y) = const. At this order, the wall 
would oscillate back and forth but would not move on 
average. Drift requires a nonzero viscosity: as the wall 
coordinate X increases on average, the loss of Zeeman 
energy must be accounted for through viscous friction. 



D. Very high field: -ff > /fo 



l^ql > 'w/2, making the steady state unreachable. As 
pointed out above, in permalloy strips with a width be- 
low 1 fim the viscous force is small in comparison with 
the gyrotropic one. As a result, the equilibrium of a vor- 
tex in the transverse direction is set mostly by the bal- 
ance of the transverse components of the gyrotropic force 
GV and the restoring force — fcl^q. The critical point is 
reached when Y^, 



w/2: 
GVc = kw/2. 



(15) 



With the aid of Eq. © we obtain the critical velocity 

(16) 



27r^ 



and the critical field 



-ffc = TxxVc/Q = kwTxx/{2QG). 



(17) 



For permalloy, 7 — 2.21 x 10^ mA^^s^^ and M — 
8.6 X 10'^ m~^ Ajii Taking the thickness of t = 20 nm we 
obtain Vc = 106 m/s. This is not too far from the critical 
velocity of 80 m/s observed by Beach et alA 

Equation ([16]) shows that the critical velocity should 
grow linearly with the film thickness t. It is easy to see 
that this result is valid beyond the crude model of a vor- 
tex wall adopted in this calculation. The two forces bal- 
ancing each other scale differently with t. While 
the gyrotropic force is linear in t, the restoring force 
comes from the magnetostatic energy, which represents 
Coulomb-like interaction of charges with density 0{t), 
hence (the dipolar part of) the restoring force is quadratic 
in t. That gives 14 oc 



C. High field: H > He- General remarks 

Numerical simulations indicate that, after the original 
vortex with a core polarization p is expelled from the 



We first demonstrate that at a very high field the ve- 
locity is again proportional to the field and calculate the 
high-field mobility. The new field scale Hq is set by the 
requirement that the restoring force —kY be negligible 
in comparison with the Zeeman force QH. The charac- 
teristic field is 



Ho = kw/{2Q) = HcG/Txx > i?c 



(18) 



When H ^ Hq, the dynamics is dominated by the Zee- 
man and gyrotropic forces, so that the vortex moves 
along an equipotential line Y = ^72 + const, or X — 2Y. 

As a result of the drift with a velocity Vd, the Zeeman 
energy goes down on average at the rate QHVd- It is 
dissipated through heat generated at the rate 



rv = 2 Or 



The transverse velocity of the vortex core reflects the 
balance between the longitudinal components of the gy- 
rotropic and Zeeman forces: Y « QH/G. We thus find 
the drift velocity 

Vd = ^{Tyy + ^Txx+^Txy)=^-^^^. (19) 

In the last transition we have used the relation between 
the coefficients of the viscosity tensor specific to this 
model 

The high-field (HF) mobility is suppressed in com- 
parison to the low-field (LF) one (fT4)) : 



Mhf 



XX 

G2 



< 1. 



(20) 



In the experiment of Beach et al* /iKp/z^LF ~ 0.1, while 
the theoretical result is 3{Txx /G)^ w 0.05, i.e. twice as 
small. 
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E. High field: H > He. Details 

To find the drift velocity of the vortex at fields above 
the vortex expulsion field H^^ we determine the total X 
displacement of the vortex over a full cycle of motion 
from the top of the strip to the bottom and back again. 
The crossing time will be slightly different on the upward 
and downward trips due to the asymmetry of the vortex 
wall and the Y component of the Zeeman force. 

Solving Eq. pT|) with polarization p = ±1 gives us the 
crossing times and displacements AT+ and AX+ (top to 
bottom) and AT_ and AX- (bottom to top): 



AX± 
AT± 



QHAT± G T TxY 



det r + G2 



fcr 



In 



XX 



i XX 

1 + HJH ± g 
1 - Hc/H ± g 



The drift velocity is 
Vd 



AT+ + AT_ 



(21) 



(22) 



(l + detr/G2)ln 



(l + ge/g) 



The resulting curve is shown for several strip widths in 
Fig. [21 Note that the critical field is not exactly He and 
actually changes slightly with the width. This is because 
the equilibrium points for both up- and downwardly po- 
larized vortices must be expelled from the strip for the 
character of the motion to change. By Eq. (|13p a down- 
wardly polarized vortex requires a slightly higher field 
to expel than an upwardly polarized one. An expansion 
of Eq. in powers of 1 /H yields the high- field result 



IV. DISCUSSION 

We have explored the dynamics of a vortex domain 
wall in a magnetic strip of a submicron width. We have 
applied the method of collective coordinates^ to the case 
when the wall has two soft modes related to the mo- 
tion of the vortex core. A simplified model of the vortex 
domain wall described in this paper yields solvable equa- 
tions of motion. The calculated mobility of the wall in the 
steady-state viscous regime at low fields agrees well with 
the value measured by Beach et alA The steady motion 
breaks down when the equilibrium position of the vortex 
moves beyond the edge of the strip. The critical velocity 



(1161) depends just on the magnetization length and the 
sample thickness; its calculated value agrees reasonably 
well with the data of Beach et alA The dynamics above 
the breakdown changes the character from overdamped 
to underdamped: the ratio of the viscous and gyrotropic 
forces acting on the wall Txx/G = 0.13 in their exper- 
iment. In this regime the velocity sharply declines at 
first but later starts to rise again as the field strength in- 
creases. The high-field mobility is reduced in comparison 
with the low-field value by the factor BF^^jf/G^ = 0.05; 
the observed reduction is not as strong: /xhf/mlf « 0.14 
In addition to simplifying the geometry (but not the 
topology) of the domain wall, we have made other as- 
sumptions that require further checking. First, we have 
assumed that any vortex absorbed by the edge is imme- 
diately reemitted. At fields lower than that required for 
emission to occur, the wall may simply stay transverse 
and continue to move in a viscous fashion. At higher 
fields, there may be short delays between absorption and 
reemission during which the motion of the wall is again 
viscous; the higher mobility of a transverse wall would 
tend to increase the drift velocity. 

Second, just as the appearance of y as a new degree 
of freedom completely changes the character of the wall 
dynamics above the critical field iJ^, at still higher fields 
additional modes of the wall may become important. The 
number and dynamical characteristics of soft modes may 
also change discontinuously as additional vortices or an- 
tivortices are created and annihilated in the bulk of the 
strip. We have observed the creation and subsequent an- 
nihilation of a vortex-antivortex pair near the original 
vortex of the wall. Like the process described by Van 
Waeyenberge et al.r^ the pair creation mediates the fiip- 
ping of the polarization of the wall vortex and results in 
the reversal of the gyrotropic force. Thus the dynamics is 
similar to that described in this paper: the vortex moves 
back and forth, while the domain wall slowly drifts along 
the strip. A possible way to detect this new regime is 
to measure the frequency of longitudinal oscillations: be- 
cause the vortex does not reach the edge, the frequency 
should be higher than 2w = 2^H expected when the vor- 
tex moves from edge to edgei^ We shall describe the onset 
of this type of motion more fully in future work. 
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